Patterning of polar active filaments on a tense cylindrical membrane 
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We study the dynamics and patterning of polar contractile filaments on the surface of a cylindrical 
cell using active hydrodynamic equations that incorporate couplings between curvature and filament 
orientation. Cables and rings spontaneously emerge as steady state configurations on the cylinder, 
and can be stationary or moving, helical or segments moving along helical trajectories. Contractility 
induces coalescence of proximal rings. We observe phase transitions in the steady state patterns 
upon changing cell diameter and make several testable predictions. Our results are relevant to the 
dynamics and patterning of a variety of active biopolymers in cylindrical cells. 



Cytokinesis, a common mechanism by which cells di- 
vide, involves the regulated assembly and constriction of 
a contractile actomyosin ring. A well studied model sys- 
tem of the dynamics of assembly of the contractile ring is 
fission yeast [1] . The contractile ring appears to assemble 
on the inner cell surface via heterogeneous nucleation of 
nodes and growth in the form of actin cables [1]. Recent 
experiments on cylindrical fission yeast and its mutants 
have shown that the stability and dynamics of the con- 
tractile ring is greatly influenced by cell geometry [2] . In 
this paper, we account for these diflPerent pheno types ex- 
hibited by the actomyosin filaments using a single model 
involving the interplay between geometry and active con- 
tractile mechanics. Indeed this interplay between geome- 
try and cytoskeletal mechanics has been well appreciated 
in the context of mitotic division in epithelial cells [3]. 

We present an extensive analysis of the steady state 
patterns of active polar filaments (e.g., actomyosin) on 
the (inner) surface of a nondeformable cylindrical cell, us- 
ing the equations of active hydrodynamics [4] generalized 
to include the coupling of geometry to filament orienta- 
tion [5]. Such nondeformable cylindrical cell surfaces are 
realized in rod-shaped bacteria and fission yeast cells, due 
to their strong coupling to a rigid cell wall [1]. While our 
study complements the approach based on microscopic 
agent-based simulations [6, 7], it has significant points 
of difference, notably the influence of cell geometry and 
size. Our main results are : (a) Rings and cables appear 
as generic steady state patterns, either from a spinodal 
instability of the homogenous conflguration or from a nu- 
cleation and growth of an actin node. Rings/cables may 
be either stationary or mobile, and the ring/cable width 
is set by a Peclet length Lp, the ratio of the fllament diffu- 
sion to a curvature-dependent active filament advection. 
(b) Phase transitions between the different steady state 
configurations can be achieved by changing cell diameter 
(Fig. 1). (c) Due to the intrinsic anisotropy in curvature 



of the cylindrical cell, the active filaments can exhibit he- 
lical shapes or tilted segments which move along helical 
trajectories, (d) Proximal rings merge to form one ring, 
the dynamics of ring merger is controlled by actomyosin 
contractility. 

We describe the dynamics of active polar filaments in 
a 2-dimensional cylindrical surface 5, in terms of a local 
concentration c(r, t) and polarization n(r, t) = ine^Uz)^ 
representing the polar orientation (equivalently, the ve- 
locity of the filaments relative to the medium). Assuming 
that momentum gets dissipated predominantly by fric- 
tion at the cortex, the hydrodynamic velocity is obtained 
from local force balance, — Fv = V-a where a = —Wcww 
is active contractile stress with 1^ < [8]. Conservation 
of total number of actin filaments in the cortex implies, 
dtc = —V • J, where, V = [R~^de^dz) is the differential 
operator on the surface of cylinder of radius R. The fil- 
ament current Ji = vocui + Aijkii^jkcni — DViC^ where 
the first term is an active advection (e.g., motor-driven 
or (de) polymerization) with speed vq^ the second is an 
anisotropic advection due to coupling between n and the 
curvature tensor /^^j, and the third is a diffusive current 
with coeflScient D. We will implicitly assume that the 
filaments are short compared to the cell size. 

The dynamics of the polarization n is obtained by gen- 
eralizing the Toner-Tu equations [9] to include symmetry 
allowed couplings to cell geometry, 

^ + A (n • V) n = V^n + K2V (V • n) 

-hCVc -h (ce - P\n\^) n -h 7 n 

(1) 

where the laplacian is defined on the surface of the cylin- 
der, and the last term implies tensor contraction. 

The 4^^ rank tensor parameters K^^^ and 7^-^^ are the 
most general phenomenological coupling between curva- 
ture tZij and polar order n to lowest order. While the 
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coupling parameter 7^-^^ can be present even in an equi- 
librium membrane with a tilt field, the parameter A^^^ is 
purely active in origin. Using symmetry considerations 
and the form of the curvature tensor hz = ( ^(^^ ) foi" 
a cylinder, we find there are two independent non-zero 
contributions from A^^^ and 7^-^^ each, denoted as A^ and 
ji with i = O^z^ henceforth. Note that A^ trivially renor- 
malize the advection in the c equation, giving rise to 
an anisotropic advection Vi = ^ ^i/R- Similarly, 7^ 
renormalize the coefficient of the linear term a in Eq. 
1 to ai = a -\- ji/R — SioKi/R'^ where 6ij is the Kro- 
necker delta. In what follows, we let 7^ to have positive 
and negative values while we restrict K12 to the positive 
domain. 

The parameter a is a linear function of the mean fil- 
ament concentration cq and governs the transition from 
orient at ionally disordered to polar ordered phase, while 
C > describes the tendency of contractile filaments to 
reorient towards the gradient of concentration [10]. Due 
to activity of the filaments they move relative to the sol- 
vent in the direction of n, thus reorientation is accompa- 
nied by movement of filaments towards each other and 
can, for instance, be enhanced by increasing motor ac- 
tivity. We explore the steady state patterns of active 
filaments on a cylindrical cell as a function of (i) cq, (ii) 
C and (iii) cell size R. Values of parameters used cor- 
respond to those in fission yeast (5. Pombe)^ and are 
tabulated in [11]. 

We first study the spinodal instability of the homoge- 
neous configuration with concentration cq, which is sta- 
ble when C = 0. When ai{co^R) < 0, i.e., when the 
filament concentration is low enough, the stable homoge- 
neous phase has no orientational order, (n) =0. We ask 
whether the homogeneous- isotropic phase (c = cq, n = 0) 
continues to be stable as we increase Using linear 
stability analysis we find that for ( > Cc({<^i}; {^i}) 
this phase become unstable with a maximally unstable 
wavevector along either the or the z axis [11], this 
gives rise to the phase boundaries depicted in [11]. The 
nature of the final steady state configurations character- 
izing the phases, however, cannot be obtained from a lin- 
ear analysis. To ascertain this we numerically integrate 
[11] the dynamical equations for c and n using an im- 
plicit and alternating direction method [12]. The finite 
wavevector instability seen in the linear analysis shows 
up as a density clumping along a specific direction. For 
parameter ranges which result in only one maxima in the 
dispersion surface, the direction of the density clumping 
is set by the direction of the fastest growing wavevec- 
tor q^rn [II], thus when is along we obtain rings^ 
while when q^ is along ^, we get cables. For the pa- 
rameter range where the dispersion surface shows two 
maxima (along both and g^), the final configurations 
correspond to asters/nodes [11]. 

When the filament concentration is high enough so that 
one or both of the a's are positive {ai{co^R) > 0), the 
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FIG. 1. (color online) Phase diagram in (^ — R ai high filament 
concentration and net mean orientation, obtained from linear 
stability analysis, for (a) ae > and ae > otz and (b) cxz > 
and az > ae- The steady state patterns corresponding 
to these phases are depicted in Fig. 2. Uniform phase refers 
to homogeneous, oriented phase along (a) and (b) z. (c) 
Scaling of the inverse of the fastest growing wave vector qm 
with C and R (inset). All quantities are in scaled units [11]. 



homogeneous solution c = cq has polar order, |(n)| 7^ 
(Fig.l). The spontaneous polarization is determined by 
the largest a and can be along either the azimuthal (^- 
axis) or the axial (z-axis) direction. Again we use linear 
stability analysis to check the stability of the uniform so- 
lution upon increasing (. Since the nature of instability 
is the same whether ao > or vice versa, we discuss, 
without any loss of generality, the former case where we 
perturb the system around no = \/ a^j ^. The system 
tends to form spatial structures only beyond a thresh- 
old contractility (" given below. As before, the orienta- 
tion of the spatial structure is set by the fastest grow- 
ing wavevector q^ of the dispersion surface : (a) when 
Cc2 > C > Cci = ^^coV.'''\ the instability is along z, 
while (b) when Cci > C > Cc2 = the instability is 

along d. Fig. 1(c) shows the dependence of the magnitude 
of the fastest growing wavevector on R and the exact 
expressions are displayed in [11]. 

Once again, linear stability analysis tells us nothing 
about the final steady state configuration, which we have 
to obtain by numerical integration of the dynamical equa- 
tions. We confirm that the linear instabilities in the para- 
mater ranges (a) and (b), give rise to periodically sepa- 
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rated stationary rings and moving cables, respectively 
(Fig. 2). However, when ( > max(Cc2, Cci), the disper- 
sion surface has two maxima, one along z and the other 
along [11], and a numerical solution of the dynami- 
cal equations, shows that the steady state configuration 
is an array of inward-pointing asters/nodes. The final 
phase diagram in is shown in Fig. 1(a), where the 
phase boundaries are obtained from linear stability anal- 
ysis, while the nature of the final steady state configu- 
rations are obtained numerically. The time sequence of 
configurations starting from the homogeneous, oriented 
phase and leading to the formation of rings is displayed 
in [11]. The phase diagram for the case when > and 
<^z > (i-e., when the spontaneous polarization is along 
z) is shown in Fig. 1(b). 

From the numerics we find that the magnitude of 
corresponds to the periodicity d of the rings, leading to 
two testable predictions (Fig. 1(c)) : (i) the distance be- 
tween rings d increases monotonically with R and satu- 
rates to a constant which depends on ( [11] and (ii) d 
decreases with contractility as 1/vT [H]- 




FIG. 2. (color online) Schematic of steady state patterns 
characterizing the phases in Fig.l : (a) stationary rings (b) 
stationary cables, (c) moving rings, (d) moving cables (e) 
asters/nodes. Small arrows (black) show texture and thick ar- 
row (blue) denotes direction of movement of the rings/cables. 

The width w of the rings and cables (and aster size) 
is set by the ratio D/vq^zi the curvature-dependent ac- 
tive Peclet length, which balances diffusion with advec- 
tion. To obtain the leading contribution to the ring 
velocity, we use the ansatz for an axisymmetric mov- 
ing ring c = c{z — ZQ{t)) together with the assumption 
that the other variables, such as w and the polarization- 
profile are fast] integration over the ring area then gives 
io(^) = VziC^R)- It is reasonable to expect > 0; 
this leads to another prediction of the model, namely the 
velocity of the moving ring is higher for narrower cylin- 
ders and approaches vq as when R is large. Similar 
arguments suggest that the velocity of moving cables is 
ve iC,R). 

We have implicitly assumed that K12 > 0, thus ob- 
taining parallel filament orientations within the ring. 
However the steady states also admit anti-parallel fila- 
ment orientations if K12 < 0; one then needs to augment 
Eq.(l) by a symmetry allowed 4th-order spatial deriva- 
tive for stability. 



Thus far we have studied the spinodal instability of the 
uniform phase to coherent structures. In the context of 
fission yeast however, actin is nucleated at the cell surface 
by the actin nucleator formin [1], which first forms a 
node, then grows into a ring or cable. We study this 
'nucleation and growth' driven transition as a function of 
Co, C R. This is easiest to study in the limit A ^ 0, 
when the the dynamics of n can be obtained from an 
"energy" functional (strictly lyapunov functional), 

E[c,n] = J (i^i+i^2)(V-n)'+i^i (V x n)'+Cc (V • n) , 

(2) 

together with a local constraint on the magnitude, n-n = 
tiq = ao^z/P' We now compare the "energies" of an 
inward- aster {Ea)^ a 6>-segment {Eq) and a z-segment 
{Ez) of fixed area A (where the segment widths are de- 
termined by the balance of the current J). We make 
the simplifying assumption that the concentration within 
these regions is uniform and so we are comparing config- 
urations with the same mean concentration cq. In one 
constant approximation, Ki = K2 = 



Ea = irKln 
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where ^ and Cc are the core size and energy, respectively. 
Note that Ea is the same as in a planar geometry, since 
the gaussian curvature of a cylinder is zero [13]. The 
"energy" branches and the phase diagram are shown in 
Fig. (3). 

Thus far, the steady state configurations correspond to 
rings (cables) along the {z axis) alone, no other orienta- 
tion is permitted. However in principle, the parameters 
accompanying the nonlocal terms in (1), namely Ki^K2 
and C, can also have anisotropics consistent with cylindri- 
cal symmetry. This could arise, for instance, from having 
"easy-directions" [14] for spatial variations of n and c. 

Consider a tilted segment at an angle to the cylinder 
axis {i/j = cos~^(m|| • z), where m|| is the unit vector 
along the segment boundary. Fig. (4b)), having a uniform 
width w{iIj) and length /(?/;) such that the area A = wl 
is fixed. For simplicity, we assume that within the tilt- 
segment, the filament concentration c = Cq is uniform. 
The texture within the tilt-segment is shown in Fig. (4b), 
the orientation of n changes by tt on a length scale w{iIj). 
The width of the tilt-segment w{iIj) is set by a balance of 
the net current J • , where is a unit vector normal 
to the boundary. In the equal constants approximation, 
Ki = K2 = the "energy" (2) of this tilt-segment is 
given by. 
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FIG. 3. (color online) Phase diagram in ^co versus (a) R 
(with A constant) and (b) A (with R constant) showing asters 
(nodes), ^-segment and 2;- segment. Insets show the "energy" 
branches (units of K) of the 3 configurations, Eq.(3b). Note 
that the form of Eg in Eq.(3b) is valid for a ^-segment which 
does not encircle the cylinder, for which the segment width is 
set by D/vz. For smaller values of R (below the kink in (a)), 
the ^-segment forms a ring, whose width is set by A/27vR. 
The orientation of n in these configurations are shown in (c). 




FIG. 4. (a) Orientation of tilt-segment schematically rep- 
resented in (b), with (short) arrows showing n within the 
segment. The steady state tilt plotted as a function of 
shows a continuous transition from a ^-segment {ip = 
to a tilt-segment (see text) at — 1. The symmetries of 
a cylinder allow for 4 stable solutions, the two shown above 
have a net orientation of n along '0, while two more have a 
net orientation of n along 



and the steady state orientation of the tilt-segment is 
obtained by setting f^|^ = 0- To prove our point, 
it suffices to look at the anisotropy of K{tlj) alone, 
which owing to cylindrical symmetry can be written as 
K{iIj) = km cos2mV^, of which we take only the first 
two modes. The phase diagram Fig. (4), shows a con- 
tinuous transition from a ^-segment with ip = ±7r/2, to 



a tilt-segment with = cos ^ y ^ ~ ^i/4A^2 

is varied. Because the net polarization of these active 
tilt-segments is along the short tilt-segments move on 
helical tracks, while the longer tilted segments appear as 
moving helices (see Refs.[15, 16] for propagating helices 
of FtsZ in B.subtilis). The analysis is relevant to recent 
experiments [17] on the movement of MreB filaments on 
helical tracks, whose strong anchoring to the cell wall 
growth machinery provides the "easy-directions" on the 
surface of these rod-shaped bacteria. 

We end with a discussion on the fate of multiple rings, 
both moving and stationary. Well separated rings mov- 
ing along the H-z-axis, result in periodic oscillations in 
time (with a time period r = L/v^, which gets smaller 
for narrower cylinders), when we impose no flux bound- 
ary conditions along z or when the cylinder is capped on 
either end by hemispheres. On the other hand, station- 
ary rings, proximal to each other, merge to eventually 
form a single ring whose width is the same as original 
rings [18]. Our analytical and numerical studies suggest 
that such ring merger occurs whenever there is a nonzero 
overlap of the concentration profiles of the rings; the dy- 
namics of merger shows an initial slow rate governed by 
K which crosses over to a faster rate given by the con- 
tractile parameters, y/Cvz (which gets larger for narrower 
cylinders). 

In summary, we have presented a detailed analysis of 
the effect of curvature-orientation coupling on the pat- 
terning of active filaments in a cylindrical cell and pre- 
dict a rich variety of steady state patterns which include 
asters (nodes), stationary and moving rings and cables. 
Our work is valid in those situations where we can ignore 
the deformation of the cell membrane, which can itself 
contribute to additional instabilities [19]. The coarse- 
grained symmetry based approach presented here is com- 
plementary to the more microscopic agent-based simula- 
tions, and has the advantage of being applicable to a 
variety of cellular contexts such as (i) actin filaments in 
fission yeast cells [1], (ii) reconstituted cytoskeletal ele- 
ments on cylindrical liposomes [18], (iii) FtsZ filaments in 
bacteria (and reconstituted in yeast [20]) and (iv) MreB 
filaments in rod shaped bacteria [17]. A more detailed 
comparison to experiments will be taken up later. Our 
work sets the stage to study the dynamics of quenches 
from one phase to another, which is appropriate to the 
study of the growth of cables or rings from the nucleation 
of nodes of acto-myosin filaments in fission yeast [21, 22]. 

We thank M. Balasubramanian and his group for dis- 
cussions and collaborations. MR thanks HFSP and NSG 
thanks the support of the Minerva Foundation (grant no. 
710589). 
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We provide a table containing the values of parameters used in the analysis and a description of 
the numerical scheme. We also provide some details of the linear stability analysis with analytic 
expressions and dispersion plots in various parameter regimes. We have also included a time sequence 
of configurations starting from the unstable homogeneous phase to the formation of rings. 



UNITS AND PARAMETER VALUES 
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We have made the dynamical equations dimensionless 
by choosing appropriate units for length, time and order 
parameter (polarization, n), 

1. Unit of length, ^ = f • Taking R = 2 /im for fission 
yeast [1], this translates to / = 0.4/im. All length 
scales are measured in units of /. 



2. Unit of time, r = Taking vq = OAjj.ms~^ for 
fission yeast [2], this translates to r = 1 s. All time 
scales are measured in units of r. 



3. Unit of magnitude of polarization, n = y^a//3. The 
values of a and /3 have been chosen to be large (= 
10) and equal so as to ensure that this magnitude 
is 1 almost everywhere. 

To obtain the values of the parameters in real units, 
we multiply the number used in the simulation by its 
appropriate dimension in units of r and n (see Table 
below) . 



Parameters 


Physical values 


Scaled values 


R 


2 /im [1] 


5(0 


L 


15 /im [1] 


37.5(/) 


D 


0.1 - l/im^s-i [3] 


0.625 - 6.25(Z) 


vo 


0.4 /ims-^ [2] 


1 (l/r) 


Ki 


0.16 fim^s~^ 


liP/r) 


K2 


fim^s~^ 




a 


10 


lO(r-i) 


li 


-20 - 10/ims-i 


-50 - 25{1/t) 


Ai 


-0.8 - 1.2/im2s-i 


-5 - 7.5(ZVr) 



Note that in the above Table, Ki has been chosen to 
be of same order as D. For simplicity, we have taken 
K2 = 0. Values of 7^ and A^ are chosen such that ^ and 
^ are of same order as a and vq. 



NUMERICAL SCHEME 

We numerically integrate the dynamical equations us- 
ing an implicit scheme with alternating direction method 
[4] with timestep At = 0.01 and gridsize Az = 0.375 and 



DISPERSION RELATIONS AND FASTEST 
GROWING WAVE VECTORS 

Low mean concentration of polar filaments 

Here we show a plot of the dispersion surface and ana- 
lytic expressions for the fastest growing wave vector when 
mean filament density is low ai < 0, i.e., when the ho- 
mogeneous phase is orient at ionally disordered. As men- 
tioned in the main text, the maxima of the dispersion 
surface is along qo or as illustrated in Fig. SI (a). 



1. When 
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the fastest 



= Ci > C> C2 = 

growing wave vector is along with a rate given 
by, 



{az-D_qlf ^Ac^VzCql (1) 



where = D + and D_ = D-K1-K2. 

The magnitude of the most unstable wavevector is 
given by the maxima ofn{qz) in (1), and is plotted 
as a function of R and in Fig. Sl(b). This takes 
a simple form when we set D_ = 0, 



qm = 



2D^ 



(2) 



It is instructive to look at the asymptotic forms of 
this equation. 
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(c) qr, 



§§^VC, for large C- 



The eigenvector corresponding to this instability 
has a large component along c, implying that this 
fastest growing instability corresponds to concen- 
tration clumping along direction {rings). This 
is verified by explicit numerical integration (see. 
Fig. S2(c) for a schematic). 
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2. When 



Dae 



, the fastest 
The expres- 



= C2 > C > Ci = 

growing wavevector is along 6 axis, 
sions for the dispersion relations and the magnitude 
of the fastest growing wavevector in this regime is 
simply obtained by interchanging z and indices of 
the parameters in (1) and (2). This instability cor- 
responds to concentration clumping along z direc- 
tion (cables). This is verified by explicit numerical 
integration (see, Fig. S2(d) for a schematic). 

3. When ( > max(Ci,C2), the dispersion surface has 
two maxima displayed in Fig. SI (a). A numerical 
solution of the dynamical equations shows that the 
steady state in this case consists of nodes (schemat- 
ically shown in Fig. S2(e)). 

The phase diagrams in the space of ((", R) appear in 
Fig. S2(a)-(b) (phase boundaries computed from linear 
stability analysis) and the corresponding steady state 
configurations are displayed in Fig. S2(c)-(e). 



High mean concentration of polar filaments 

Here we show a plot of the dispersion surface and ana- 
lytic expressions for the fastest growing wave vector when 
mean filament density is high and uniform, i.e., either ao 
or az is positive, corresponding to orient ational order 
along or z. As mentioned in the text, without any loss 
of generality, we consider the case when ag is positive 
and greater than In this case, the ordering is along 
the ^-axis. 



As before, the maxima of the dispersion surface is along 
qe or g^, as illustrated in Fig. S3(a),(b), 

1. When = Ce2 > C > Cci = ^^f^, the fastest 
growing wavevector is along z^ with a growth rate 
is given by. 



= 



az-ae- D^ql 1 



(3) 



where D± have been defined before. The magni- 
tude of the most unstable wavevector in the limit 
D_ = 0, is given by 




which has the following asymptotic forms. 
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VT for large C (see Fig. 1(c) in main text). 

2. When Cci > C > Cc2, the most unstable wavevector 
is parallel to 6>-axis. The dispersion relation for the 
unstable mode and the magnitude of the fastest 
growing wavevector are given by. 



-2ae + D+ql 



-2ae + D-ql + i 




A) + AcQVeCqe 



(5) 



The magnitude of the fastest growing wavevector 
in the limit D_ = is. 



2D. 



(^J{ve - A)2 - 4coCt^,)' - l&Dlal 



(6) 



The final steady states corresponding to (a)-(b) are 
obtained from a numerical solution of the dynamical 
equations and correspond to stationary rings and 
moving cables. When > max((^ci5 Cc2)5 the numerical 
solution shows that asters/nodes are the final steady 
state configurations. The phase diagrams appear in 
Fig. l(a)-(b) in the main text. 



TIME SEQUENCES OF CONFIGURATIONS 
FROM NUMERICAL INTEGRATION 



In Fig. S4 we show snapshots from the numerical inte- 
gration of the dynamical equations when filament density 
is high, such that ae > max(a;2,0) and Cc2 > C > Cci- 
The initial homogeneous, orient at ionally ordered config- 
uration along ^, is immediately unstable to the formation 
of rings. This spinodal instability happens exponentially 
fast. The subsequent dynamics is a slow coarsening of 
the rings, which involves ring diffusion and coalescence, 
resulting in a reduction in the number of rings over time. 
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FIG. S2. (color online) (a)-(b) Phase diagrams in ^ — ^ at low filament concentration, a^ < 0, corresponding to zero mean 
orientation (n) = 0. The final steady state configurations obtained from a numerical solution of the dynamical equations for c 
and n, correspond to (c) rings, (d) cables and (e) nodes, (a) When > C2, the homogeneous phase is unstable to rings and 
nodes as ( increases, (b) When C2 > C15 the homogeneous phase is unstable to cables, at large R. As R is lowered, the phase 
boundaries cross and the sequence of transitions resembles (a). 
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FIG. S3, (color online) Dispersion surface when mean density of filaments is high and the homogeneous phase is orientationally 
ordered along direction. The maxima of dispersion surface is either along (a) 2;-axis or (b) both z and ^-axes. 
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FIG. S4. (color online) Time evolution of configurations starting from an (a) initial homogeneous, orientationally ordered 
configuration along ^, which shows an (b) early time spinodal instability toward ring formation, followed by a (c) slow coarsening 
regime. Colorbar denotes the concentration of actin filaments, in units of mean concentration cq. 



